Experimental observation of spectral gap in microwave n-disk systems 
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We present experimental studies of the symmetry reduced three-disk and five-disk systems using 
a microwave setup. By extracting the complex resonances from the signal by means of the method 
of the harmonic inversion we can access the distribution of the imaginary parts or the width dis- 
tribution. While increasing the opening of the systems, a spectral gap is observed for thick as well 
as for thin repellers and for the latter case it is compared with the known topological pressure 
bounds. Furthermore the maxima of the distributions are found to coincide for a large range of the 
distance to radius parameter with half of the classical escape rate. These results confirm theoretical 
predictions based on rigoros mathematical analysis in the case of the spectral gap and on numerical 
experiments in the case of the maxima of the distributions. 
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In semiclassical physics one is interested in the re- 
lation between quantities from classical mechanics and 
from quantum mechanics asymptotically for small ef- 
fective Planck constant of the system. Two manifesta- 
tions of the classical quantum correspondence for closed 
systems are the Weyl law [l[ and the Gutzwiller 0, Q 
trace formula. The Weyl law describes asymptotically 
the number of quantum mechanical energy levels in terms 
of the classical phase space dimension. The trace formula 
describes the fluctuations of the density of states around 
the Weyl law by means of classical periodic orbits and 
their stability [3|. 

For open systems there are also results connecting clas- 
sical and quantum quantities [U, H| ■ Their derivation is 
however much more delicate because the classical energy 
shells in phase space become non compact and the real 
valued eigenvalues of the Hamiltonian are replaced by 
complex resonances For open systems the imag- 

inary parts of resonances are always negative and they 
correspond to the rate of decay of the corresponding un- 
stable states. 

A celebrated example of a semiclassical formula is the 
fractal Weyl law for open chaotic systems with fractal 
repeller. It gives the asymptotic behavior of the counting 
function, that is the number of resonances up to some 
energy, for a fixed bound of the imaginary part. There 
exists numerous mathematical [j| [l(| , numerical (TTl - fl4j 
and experimental [IB-[l3] studies during the last years. 
Another line of research focuses on the distribution of 
the imaginary parts in the semiclassical limit [14j. 

A paradigmatic open system with a classical fractal 
repeller is the n-disk system (see Fig. [T]). It has been 
introduced in the late 80s by Ikawa in mathematics [l8T | 




FIG. 1. (Color online) Left: Sketch of three-disk system, 
where one fundamental domain is shaded. Right: Photograph 
of the experimental cavity without top plate including the disk 
inset and the absorber. 



and by Gaspard, Rice [19|, l20j and Cvitanovic, Eckhardt 
[2lj j in physics. A n-disk system describes the scattering 
of a classical or quantum mechanical particle in a two- 
dimensional plane on n hard disks whose centers form 
a regular polygon. The distance between the centers is 
historically denoted by R and the disk radius by a. Up 
to scaling the geometry is completely defined by the ratio 
R/a (see Fig. [J). 

The quantum system is described by the Helmholtz 
equation 
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with Dirichlet boundary conditions at the disks bound- 
ary. As the system is open the quantum resonances 
k n = Re fc n + ilm fc„ are complex valued. For the three- 
disk system Gaspard and Rice [Hf gave an explicit ex- 
pression for the scattering matrix in terms of Bcssel and 
Hankel functions which allows to calculate the quantum 
resonances numerically. 

Classically one studies the trajectories of particles 
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which perform hard wall reflections on the disks. From 
periodic trajectories a wide range of classical quantities 
like the classical escape rate, the fractal dimension of the 
repeller and the topological pressure can be calculated 
using the classical Ruelle (-function 20] , 
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where the product runs over the primitive periodic or- 
bits, T p are the corresponding period lengths and A p the 
stabilities. The topological pressure P(/3) is then defined 
as the smallest real pole of C/s{ z )- An effective way of 
its calculation is the cycle expansion, introduced by Cvi- 
tanovic and Eckhardt [2l| [23[ . The classical escape rate 
is given by 7 = — P(l) and the reduced Hausdorff di- 
mension dn of the fractal repeller by the Bowen pressure 
formula P{d H ) = H). 

In 1988 and 1989 Ikawa and Gaspard-Rice indepen- 
dently derived a connection between the topological pres- 
sure, being a purely classical quantity, and the quantum 
mechanical spectral gap. A spectral gap in this context is 
a constant C < such that Im k n < C. Gaspard and Rice 
started from Gutzwillcr's trace formula and used semi- 
classical zeta functions to conclude that Imfc ra < P(l/2). 
They also confirmed this estimate numerically [l9| . Later 
this estimates for the spectral gap has been also obtained 
mathematically for more general semiclassical systems 
25]. It is, however, known that this bound is in gen- 
eral not optimal as it does not take into account phase 
cancellation (see section 8.2 in [5j and [26|). Especially 
for weakly open systems this bound has no implication: 
As P(/3) is monotonously decreasing it is evident that for 
systems with dn > 1/2 the bound P(l/2) is positive, but 
the imaginary parts are negative anyway. Such systems 
are nowadays also called systems with "thick repeller" in 
contrast to systems with "thin repeller" where dn < 1/2 
and P(l/2) is a true constraint on the imaginary parts 
(see e.g. [H) . 

The same estimate on the spectral gap has been ob- 
tained in the 70s by Patterson and Sullivan [27|, [28j for 
scattering on hyperbolic quotients T\H 2 of infinite vol- 
ume which provide a mathematical model for chaotic 
scattering [29|]. In that case, quantum resonances (de- 
fined as poles of the scattering matrix of the surface) 
coincide with the zeros of the Selberg zeta function and 
the topological pressure can be calculated explicitly us- 
ing S, the dimension of the limit set of F: P(/3) = 6 — /3. 
Results of [13, [H[ imply, that all resonances satisfy 
Imk < 5 — 1/2. For this model the estimate is known to 
be sharp as there always exists a resonance at i{8 — 1/2) 
(corresponding to a bound state when 8 — 1/2 > 0). How- 
ever, there are no other resonances for Im k < P(l/2) — e, 
for some small e > (30j . The question of further im- 
provements for the spectral gap is an active field of math- 
ematical research with deep applications to number the- 
oryiHIl. 

Another interesting property of the Im fc-distribution 
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FIG. 2. (Color online) Three-disk system at R/a = 2.88: 
in the lowest panel resonances belonging to good (orange 
crosses) and the best reconstruction (blue triangles) are shown 
in the complex plane within a small frequency range, the up- 
per 2 panels include the measured signal (black, dashed) and 
best reconstruction (blue) in this window (vertical lines) based 
on the poles marked by the blue triangles, presented in real 
and imaginary part 



has been observed in [T3[ : The imaginary parts of numer- 
ically calculated resonances were found to concentrate at 
Imfc = — 7d/2 = P(l)/2 that is at the half of the classi- 
cal escape rate. This corresponds to the classical expecta- 
tion as the quantum mechanical probability density is the 
modulus square of the wave function [191 ]. Although no 
mathematical result supports the claim that imaginary 
parts of resonances concentrate at —j c \/2, Naud [32j has 
shown that the density of resonances for Imfc > — 7 c i/2 
is lower than the prediction from the fractal Weyl law 

Another connection between the classical escape rate 
and the quantum spectrum has been experimen tally ob- 
served in microwave studies of n-disk systems |l5l . [l6| : 
Lu et al observed that the decay of the wave- vector auto- 
correlation function for small wave vectors is related to 
the classical escape rate fl5) . 

In this Letter we will focus on the spectral gap 
and the maximum of the experimentally obtained Im fc- 
distribution and compare them to the calculated topolog- 
ical pressure and the classical escape rate, respectively. 
We would like to note that all following investigations 
were performed in the low lying fc-regime. 

In order to simplify the n-disk systems we exploit their 
D n symmetry and study the symmetry reduced systems. 
In the reduction to the fundamental domain also the two 
enclosing symmetry axes are hard wall potentials acting 
as "mirrors" (see shaded area in Fig. [T]). For the quan- 
tum mechanical system the Dirichlet boundary condition 
at the symmetry axes imply that the corresponding scat- 
tering resonances belong to the A2 representation. 

In the experiments we realized the quantum version of 
the symmetry reduced three- and five-disk system using 
a microwave cavity. The triangular resonator (Fig. [T]) 
has two metallic side walls of length 100 cm including an 
angle of 60° for the three-disk, and 36° for the five-disk 
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FIG. 3. (Color online) In the left panel the resonances for R/a — 5.5 in the complex fc-plane are shown as well as the distribution 
of the imaginary parts of k in the right panel. The shown k range corresponds to a frequency range from 2 to 24 GHz. The 
orange clouds correspond to all resonances leading to good reconstruction as well as the orange histogram. Note that many 
orange dots might overlay each other. Blue triangles and blue dashed-dotted histogram describe the set belonging to the best 
reconstruction. Black crosses are the numerically calculated poles and the solid black histogram the corresponding distribution. 
The dotted red line in the right panel is P(l/2), the red dashed line — 7ci/2. 



system, respectively, and the third side covered with ab- 
sorbers avoiding reflections at the open end. The R/a 
parameter can be changed by moving a half-disk inset 
of radius a — 19.5 cm, which is in contact with top and 
bottom plate, along the side wall. In the three-disk sys- 
tem the R/a range between 2.26 and 6.17 was technically 
accessible, in the five-disk case we were able to measure 
between 2 and 3.9. The half-disk inset was moved in steps 
of 10 mm. A short wire antenna (r=0.7mm) is inserted 
through a hole in the top plate, placed in the interior 
region, and it is sufficiently short not to touch the bot- 
tom plate. The height of the cavity h — 6 mm leads to a 
cutoff frequency of 25 GHz. Within the frequency range 
from 2 to 24 GHz only the TMo mode can propagate and 
the cavity may be considered as two-dimensional. Hence 
the equivalence between wave mechanics and quantum 
mechanics, i. e. between the time independent Helmholtz 
and Schrodinger equation, exists (for details of the ex- 
perimental setup see (l~7| . for a general introduction to 
microwave billiards see chapter 2.2 in (36[). 

A vector-network-analyzer (VNA) is used for measur- 
ing and reveals the complex S'-matrix. Under the as- 
sumption of a point-like antenna the measured reflection 
signal equals [37 [ 



S n (v) = 1 
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where Vj are the complex valued resonance positions. 
The extraction of Vj and Aj from the signal is the task 
of the data analysis: For closed systems and low fre- 
quencies the resonances are well separated and could be 
treated by a multi-Lorentz-fit. But for open systems, 
where the resonances overlap strongly, a fit would not 
converge. Therefore we applied the harmonic inversion 
(HI) [38l - l40l ] on the signal, a sophisticated nonlinear algo- 
rithm, to extract complex poles from the measured sig- 
nals. To guarantee an equal treatment we used it for all 
measurements. The HI was already used to analyze ex- 



perimental line widths in the presence of wall absorption 
and additional open channels [4(j and to study the frac- 
tal Weyl law [17]] • In the latter our way of using the HI 
was already discussed in detail, thus we state here only 
the main ideas: For non-ideal, i.e. experimental, data we 
have shown that the HI should be applied several times 
with different sets of its internal parameters, any leading 
to a set of complex poles and amplitudes. Afterwards the 
reconstruction based on these results is compared to the 
original signal. Fig. [5] shows part of a typical spectrum 
(black solid line) and the best (concerning the x 2 -error) 
individual reconstruction (blue dashed line) within the 
window indicated by the vertical lines. The correspond- 
ing resonances are marked by blue triangles in the lower 
panel, the complex plane. The orange crosses belong to 
other resonance sets, also leading to good reconstructions 
(to maintain clarity they are not shown in the upper two 
panels), called good resonances. Other sets not meeting 
the criterion are rejected for further analysis. 

For the three-disk system we checked the reliability 
of the HI by comparing the experimental resonances 
with calculations based on the algorithm of Gaspard and 
Rice (22[. However even for experiments with closed mi- 
crowave systems it is known that only the lowest lying 
resonances agree well with the theoretical predictions. 
For higher frequencies the experimental perturbations 
disturb the measured spectrum such that the measured 
resonances cannot be associated directly to the theoreti- 
cal ones, but statistical properties such as the resonance 
density persist. For an open system we expect this effect 
to be even more significant as the spectra of non hermi- 
tian Hamiltonians are known to be much more unstable 
under perturbations compared to hermitian Hamiltoni- 
ans mi. 

Fig. [3] shows again the good Hl-resonances in orange 
and the best in blue for R/a = 5.5 and in the full mea- 
sured Refc-range from 40 to 500 m _1 . Again the or- 
ange poles form "clouds" around the blue triangles - 
the elongated shape of the clouds is a consequence of 
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the non-isometric axis ranges. The black crosses indicate 
the numerically calculated resonances. The composition 
of resonance chains is typical for large R/a parameters 
[IH HH Ell- The individual resonances are not repro- 
duced by the experimental data due to inevitable reflec- 
tions at the absorbers and the perturbation by the an- 
tenna but the resonance free regions and the resonance 
density coincide. 

In the right part of Fig. [3] the corresponding Imfc dis- 
tributions P(Imfc) are shown, in solid black for the nu- 
merically calculated and in dashed-dotted blue for the 
experimental spectrum. Both distributions show up to 
be the same within the limits of error. The same was true 
for all good reconstructions passing the \ 2 criterion. All 
P(Im k) distributions shown below have hence obtained 
by superimposing the results of all good reconstructions. 
One example is shown in the right hand part of Fig. [3] 
in orange. In fact one can show that P(Im k) is robust 
with respect to errors in the reconstruction as long as the 
number of resonances entering the reconstruction is ap- 
proximately the same. For the example shown in Fig. [3] 
the number varied between 94 for the numerical data and 
117 for the individual reconstruction. 

By measuring the averaged Im k distribution for sym- 
metry reduced three- and five-disk systems for different 
R/a parameters we can now study the dependence of the 
Im k distribution on the opening of the system (Fig. 0]) . 
For every R/a value the averaged histogram of the Im k 
distribution is plotted as a shade plot. For the five-disk 
case the measured data is presented in Fig. Ufa). For 
R/a = 2 the system is completely closed however we ob- 
serve already a small gap m 0.15 m _1 caused by antenna 
and wall absorbing effects. While opening the system 
the very narrow Im k distribution first gets wider and 
the maximum of the distribution moves towards higher 
imaginary parts. From R/a 2.5 the resonance free re- 
gion starts to grow and reaches a value of « 0.5 m -1 for 
the maximal accessible opening at R/a = 3.9. Over the 
whole R/a range the value of P(l/2) stays positive, thus 
providing no lower bound on the spectral gap. The solid 
black line in the shade plot shows half the classical es- 
cape rate calculated by the cycle expansion. We show 
this curve down to R/a — 2.41, because for lower R/a 
values so called pruning sets in (for orbits up to order 4) : 
Due to those periodic orbits which are destroyed when 
the disks are too close the symbolic dynamic is not com- 
plete an ymo re -which was fundamental in the cycle ex- 
pansion [2l[ . In accordance with calculations in the high 
frequency regime [l3j our experiment in a much lower 
frequency regime shows that the maximum of the Im k 
distribution is described by — 7 c i/2. We emphasize that 
there are no free parameters to fit j c \ to the experiments. 

A more open systems is the three-disk system. Here 
the repeller becomes thin for R/a > 2.83, i.e. P(l/2) 
provides a lower bound on the gap. Again one observes 
that the gap first increases and only for high R/a val- 
ues coincides with the lower bound P(l/2) (dotted black 
line). At which R/a value exactly the gap in the ex- 
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FIG. 4. (Color online) Shade plots of the distribution of Irak 
using a color code from white to dark blue for (a) the five- 
disk experiment, (b) the three-disk experiment, and (c) the 
three-disk simulation as a function of the R/a parameter. The 
classical escape rate (precisely: — -y c i/2) is the solid black line. 
P(l/2) corresponds to the dotted line. For the five-disk case 
the P(l/2) is still positive and hence does not arise in the 
plot. 



periment opens and whether it opens before P(l/2) be- 
comes negative is less clear as in the five-disk system. 
The area of the totally closed system (R/a — 2) is too 
small to enable meaningful measurements. Since there 
are no pruned orbits until order 4 from R/a = 2.01, we 
are able to plot the calculated curve for the full measured 
range. The maximum of the Im k distribution decreases 
for R/a > 3 which might be surprising at first sight. The 
reason is that the time of flight between two scattering 
events increases linearly which will overcompensate the 
defocussing effect of a scattering event for large enough 
R/a. This leads to the decrease of —j c \/2 (solid black 
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line) , following the maximum of the distribution. 

As a consistency check Fig. @Jc) shows the shade plot 
for the numerical data of the symmetry reduced three- 
disk system. Also there the correspondence of — 7 c i/2 
(solid black line) is clearly visible only for large and less 
prominent for smaller R/a values. For large R/a the 
lower bound P(l/2) (dotted black line) coincides well 
with the numerically observed gap. The opening of it 
however is not described by P(l/2) which becomes obvi- 
ous for 2.5 < R/a < 2.83. Here P(l/2) is still positive 
but a clear gap is already visible, the same phenomenon 
as observed in the experimental data of the five-disk sys- 
tem. 

In this Letter we have demonstrated the existence of a 
spectral gap in open chaotic n-disk microwave systems. 
We could extract the resonances from the measured sig- 
nal and thus had direct access to the gap and the max- 
imum of the Im k distribution. These were compared 



with the calculated classical values for P(l/2) and 7ci/2. 
A good agreement was found for sufficiently open sys- 
tems. But we could also show that the bound P(l/2) 
does not describe the opening range of neither the exper- 
imental nor the exact quantum mechanically calculated 
data. We would like to emphasize that all investigations 
were performed in the low lying k regime thus showing a 
remarkable agreement with the semiclassical predictions. 
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